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ABSTRACT 
A  two-dimensional  dynamic  model  of  a  towed  buoy  designed  to  remain 
submerged  just  below  the  water  surface  is  developed.    The  study  includes 
the  model  development,  the  analysis  of  certain  model  parameters,  the 
development  of  a  tension  force,  and  the  range  of  stability  of  the  model  to 
initial  condition  perturbations  . 
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I.     INTRODUCTION 

The  purpose  of  this  investigation  is  to  attempt  to  mathematically 
model  and  computer-simulate  the  two-dimensional  motion  of  a  submarine- 
towed  buoy.     This  has  a  very  specific  military  application  in  submarine 
communications  systems;  the  basic  system  configuration  as  a  submarine 
communication  link  is  illustrated  in  Fig.    1.     It  is  desired  that  the  buoy 
stay  underwater  at  all  times,  and  remain  within    ±1  ft.  of  four  feet  below 
the  water  surface;  this  gives  the  buoy  an  operating  range  from  three  to 
five  feet  below  the  surface.     An  additional  system  requirement  is  that  the 
buoy  system  stay  within  its  allowable  operating  zone  while  following 
waves  up  to  sea-state  four. 

Many  investigators,   such  as  Lacey,  have  dealt  with  ship-towed  or 
airplane-towed  devices,  where  the  buoy  is  the  "bottom-following"  type. 
Pode  and  other  early  investigators  studied  the  steady-state  configuration 
of  cables  in  a  steady  uniform  stream;  and  more  recently  Schram  developed, 
by  the  method  of  characteristics,  a  two-dimensional  dynamic  cable-buoy 
model.     However,  he  treated  the  buoy  as  a  vector  force  and  did  not  go 
into  the  actual  buoy  dynamics. 

Steady-state  studies  do  not  show  the  response  of  the  buoy  and  cable 
to  perturbations  such  as  sea  state.     It  is  this  response  that  will  prove 
most  useful  in  the  study  of  buoy  and  cable  parameters.    This  study  has 
been  primarily  concerned  with  develo  ing  the  dynamic  model  of  the  buoy; 
and  the  tension  force  generated  by  the  cable  has  been  loosely  approximated. 
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The  model  is  of  a  buoy  presently  being  developed  for  the  U.  S.  Navy; 
thus,  the  dimensions  and  other  data  were  obtained  from  an  existing 
prototype. 

This  study  develops  the  two-dimensional  model  around  two  basic 
equations;  force  equals  mass  times  acceleration,  and  the  moment  equals 
the  mass  moment  of  inertia  times  the  angular  acceleration.    The  model 
block  diagram  was  built    around  these  two  equations,  several  basic 
hydrodynamic  relationships,  and  the  wind-tunnel  and  dimensional  data 
supplied  on  the  buoy.     Once  the  block  diagram  was  formulated,  it  was 
then  ready  to  be  computerized.     The  model  was  simulated  using  the    IBM 
3  60-67  at  the  Naval  Postgraduate  School,  and  DSL  (Digital  Simulation 
Language);  DSL  is  a  very  simple,  but  versatile  simulation  language 
developed  by  IBM.     Once  the  model  was  computerized,  the  lengthy  and 
time-consuming  task  of  optimizing  certain  parameters  and  finding  a  real- 
istic formulation  of  the  tension  force  followed.     The  determination  of  the 
tension  force  and  its  coupling  with  the  buoy  has  been  one  of  the  most 
difficult  problems  in  the  formulation  of  the  buoy-cable  model.     This  has 
been  a  problem  in  all  buoy-cable  models  where  the  dynamic  response  is 
of  interest;  this  is  mainly  due  to  the  fact  that  it  is  not  exactly  known  how 
the  forces  generated  by  the  buoy  and  cable  react  on  and  with  each  other. 
The  development  of  the  model  and  the  results  of  the  simulation  are 
presented . 
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II.     EQUATIONS  OF  MOTION  AND  PROBLEM  SETUP 
The  purpose  of  the  buoy  is  to  remain  four  feet  below  the  water 
surface  so  that  a  protruding  antenna  may  maintain  two-way  communica- 
tions; when  the  buoy  is  operating  in  high  seas  it  must  follow  the  sea 
surface.    The  body,  as  shown  in  Fig.  2.,  is  basically  a  fat  wing  that 
produces  dynamic  lift  according  to  its  angle  of  attack.     This  angle  of 
attack  may  be  changed  by  torque  produced  by  the  fin.     Besides  the 
dynamic  lift,   produced  by  the  body  angle  of  attack,  there  is  the  static 
lift  of  positive  buoyancy.     The  buoy  model  is  based  on  the  two  basic 
equations: 

F=ma 
M=I    V 

y 

The  first  equation  is  used  to  describe  movement  in  the  horizontal  and 
vertical  planes,  the  X  and  Z  planes  respectively.     The  second  is  used 
to  describe  the  rotation  of  the  body  about  the  Y  axis;  these  axis  are 
shown  in  Fig.  2.    The  three  hydrodynamic  forces  of  lift,  drag,  and  moment 
are  defined  as  follows: 
L=£CLpV2S 

D=KDpV2S 


2      Mp 
Along  with  the  tension  force  produced  by  the  cable  and  the  buoyancy  of 
the  body,  these  three  equations  describe  all  the  forces  acting  on  the 
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body.     Before  going  into  the  complete  buoy  block  diagram,  the  body 
coefficients  will  be  put  into  a  more  convenient  form,  and  the  fin  servo 
loop,  which  controls  the  buoy  attitude,  will  be  examined. 

A.      BODY  COEFFICIENTS 

In  both  aerodynamics  and  hydrodynamics  ,  there  are  various  dynamic 
coefficients;  these  are  determined  by  the  shape  of  the  body  involved. 
There  are  tables  of  coefficients  for  certain  common  wing  shapes.     However, 
when  an  odd  body  configuration  is  being  studied,  it  may  not  be  possible 
to  apply  the  standard  shapes  in  order  to  derive  overall  coefficients  for 
the  body.     In  this  case  the  coefficients  of  the  body  must  be  determined 
by  wind-tunnel  tests  . 

This  was  the  case  with  the  system  being  studied;  the  body  is  of 

non-standard  shape,  but  the  fin  is  the  standard  NACA  0012  shape.    The 

three  coefficients,  CT  ,  C^,  and  C,„,  were  measured  at  various  body 

L        D  M 

and  fin  angles.    The  coefficient  of  lift  and  the  coefficient  of  moment 
include  the  effect  of  the  positive  buoyancy  force  of  490  pounds.    The 
effect  of  the  antenna  drag  is  not  included  in  any  of  the  coefficients  . 
The  values  of  these  coefficients  are  shown  in  Tables  I— III  as  a  function 
of  body  angle  of  attack  and  fin  angle  of  attack.     It  was  decided  that  for 
computer  implementation,  these  coefficients  would  be  easier  to  handle 
if  they  were  in  equation  form.     This  was  done  by  applying  a  least-squares 
fit  twice  to  each  data  set;  the  resulting  equations  may  be  found  in  Sub- 
routine Tink  of  the  digital  program  of  the  model.     The  resultant  curves  and 
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FIN  ANGLE  (Deg.) 

-10     -7.5      -5.     -2.5  0.0  2.5  5. 

-12  -.054    -.039    -.027    -.016  -.003  .0064  .0211 

-9   -.001     .0110    .0253    .0353  .0477  .0589  .0756 

-6    .0585    .0647    .0776    .0901  .1020  .1139  .1306 

-3    .1102    .1214    .1351    .1461  .1593  .1723  .1876 

BODY  0    .1676    .1792    .1929    .2035  .2166  .2294  .2454 

ANGLE 

(Deg.)  3    .2315    .2432    .2573    .2677  .2807  .2943  .3095 

6    .3044    .3163    .3285    .3414  .3547  .3678  .3828 

9    .3834    .3946    .4075    .4195  .4333  .4461  .4593 

12   .4654    .4774    .4911    .5025  .5155  .5298  .5386 

COEFFICIENT  OF  LIFT 

Table  I 
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BODY 

ANGLE 

(Deg.) 


FIN  ANGLE  (Deg.) 

-10     -7.5     -5.      -2.5     0.0  2.5     5. 

-12  .0415    .0374    .0364    .0331    .0327  .0307  .0320 

-9   .0345    .0314    .0310    .0284    .0283  .0271  .0288 

-6   .0319    .0302    .0298    .0279    .0287  .0284  .0306 

-3   .0347    .0334    .0336    .0322    .0343  .0347  .0377 

0   .0403    .0404    .0415    .0404    .0430  .0444  .0475 

3   .0514    .0522    .0537    .0534    .0568  .0583  .0625 

6   .0696    .0703    .0721     .0730    .0769  .0804  .0844 

9   .0956    .0972    .0998    .1007    .1048  .1087  .1143 

12  .1306    .1326    .1354    .1377    .1428  .1475  .1522 

COEFFICIENT  OF  DRAG 


Table  II 
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FIN  ANGLE  (Deg.) 

-10     -7.5     -5.      -2.5      0.0  2.5     5. 

-6   .0981    .0908    .0834    .0772    .0707  .0639  .0554 

-3   .0794    .0723    .0651     .0593    .0528  .0452  .0370 

0   .0585    .0512    .0439    .0382    .0314  .0237  .0158 

BODY  3   .0539    .0289    .0214    .0156    .0086  .0010  -.007 

ANGLE 

(Deg.)  6   .0106    .0036   -.003    -.009    -.017  -.024  -.032 

9  -.016    -.023    -.030    -.036    -.044  -.051  -.058 

12  -.044    -.050    -.057    -.064    -.071  -.078  -.-84 

COEFFICIENT  OF  MOMENT 
Table  III 
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data  points  for  constant  body  angle  of  attack  and  varying  fin  angle  of 
attack  are  shown  in  Figs  .   3-5  . 

B.      FIN  SERVO  LOOP 

The  diagram  of  the  fin  and  the  servo-control  system  is  shown  in 
Fig.   6.    The  purpose  of  the  fin  is  to  orient  the  body  so  that  it  will  stay 
four  feet  below  the  water  surface.    The  fine  produces  a  moment  on  the 
buoy  which  causes  the  body  angle  of  attack  to  increase  or  decrease. 
This  increase  or  decrease  in  body  angle  of  attack  causes  the  body  lift 
to  increase  or  decrease  and  the  body  rises  or  descends  as  is  necessary 

The  fin  actuator  loop  is  driven  by  a  Size  23,  Kearfott  motor.     The 
loop  is  designed  to  have  a  rise  time  of  0.04  sec;  this  actuator  loop  is 
shown  in  Fig.   7.    The  servo  loop  used  in  the  buoy  model  is  considered 
to  be  a  K/s2  plant  with  tachometer  and  unity  feedback;  this  is  shown  in 
Fig.  8.    The  system  transfer  function: 

G  KA /s* 


1    +     GH  1    +  K  /s2(K  s  +  1) 


KA 


sa  +  K.K.s    +  K, 


A  t  A 

and  the  characteristic  equation: 


s2  +   K.K.s    +  K      =   0 


A  t  A 


Assuming  K  =1000 
A 


w2  =  KA  =  1000 

n    A 
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s     n  At 


24w 
K     =  n      . 

1000 
Taking  a  value  of    4=0.7: 

K     =   2(0.7)(31.62)/1000 

=    .0443     . 

With  K„    and  K     set  at  the  above  values  ,  the  rise  time  of  the  servo 
A  t 

loop  is  about  0.1  sec.     This  is  not  as  fast  as  desired,  but  by  changing 
the  two  gains  the  desired  rise  time  may  be  obtained.     Two  different 
responses  of  the  fin  servo  loop  are  shown  in  Fig.  9;  both  have  an  input 
step  of  0 .  174533  rad.    (10  degrees).    The  first  set  of  gains  is  used  in  the 
model  simulation.     This  is  because  it  is  doubtful  that  the  inputs  into  the 
servo  loop  will  require  that  the  fin  move    that  fast;  if  so,  then  the  buoy 
has  probably  reached  the  limits  of  its  surface-following  capability. 
The  inputs  into  the  fin  servo  loop  are  shown  in  Fig.  6.    The  fin 
servo  loop  that  was  simulated  above  is  composed  of  the  summer,  servo 
amplifier,  motor  and  gears,   linkage  to  the  fin,  and  the  unity  and  tacho- 
meter feedback  from  the  fin  to  the  summer.    The  other  three  inputs  will 
come  from  gages  and  electronic  sensors  in  the  body.     These  inputs  are 
the  depth  error,  the  depth  rate,  and  the  body  pitch  rate.     In  the  model 
simulation,  there  are  three  constants  (ADE,  ADR,  and  APR)  which  multiply 
these  three  inputs.    The  values  of  these  three  constants  are  discussed 
later  in  the  study. 
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III.     THE  BUOY  MODEL 
In  this  development  of  the  buoy  model,  the  effect  of  water  currents, 
the  forces  generated  by  the  heaving  water  surface,  and  the  free-surface 
effect  are  not  taken  into  account.     Currents  can  be  considered  neglible 
compared  to  the  heaving  forces  the  buoy  might  experience.    The  free- 
surface  effect  is  experienced  by  the  buoy  because  it  is  a  body  moving  at 
a  depth  comparable  to  its  length;  a  more  detailed  and  complicated  model 
would  have  to  include  this  force.     The  force  generated  by  the  heaving 
water  surface  would  not  affect  the  buoy  as  much  as  might  be  expected. 
This  is  because  the  body,  as  long  as  its  vertical  freedom  is  not  restricted, 
will  act  like  a  particle  in  the  water  and  will  move  up  and  down  with  the 
surface.    Thus,  the  heaving  forces  will  actually  help  maintain  the  buoy 
near  the  water  surface,  which  is  what  is  desired.     If  the  buoy  were  being 
towed  directly  above  its  towpoint,  it  would  have  little  or  no  vertical 
freedom;  however,  if  the  buoy  is  well  behind  the  towpoint,  then  the  cable 
will  provide  little  restriction  to  the  buoy's  vertical  freedom.     Tests  on  the 
actual  buoy  have  shown  that  with  475  ft.  of  cable  and  the  two  point  at  a 
depth  of  2  40  ft.  ,  the  buoy  can  heave  up  to  2  0  ft.  with  the  waves  without 
the  cable  decreasing  the  buoy's  vertical  freedom  of  motion  by  more  than 
10  percent.1 


1  Soulant,  H.,     Sea  Surface  Followage  by  Towed  Buoy  (Preliminary 
Study),   p.    II,     M.irch  1969. 
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The  buoy  model  will  now  be  developed;  the  buoy  diagram  is  again 
produced  with  some  specific  dimensions   (Fig.   10).     Besides  the  diagrammed 
dimensions,  the  buoy  has  the  following  additional  characteristics: 

S   =   14.25  ft2 

C  =      6  . 0  f  t 

b   =      2.375  ft 

mass  =  70  slugs 

mass  moment  of  inertia  =  975  slug -ft2 

volume  =  17.5  ft3 
The  block  diagram  of  the  dynamic  model  is  shown  in  Fig.   11.    As 
was  previously  stated,  the  buoy  model  is  based  on  the  summation  of 
forces  in  the  X  and  Z  planes  ,  and  about  the  Y  axis  .    The  summed  forces 
in  the  X  and  Z  directions  are  divided  by  the  mass  of  the  buoy,  and  the 
moments  about  the  Y  axis  are  summed  and  divided  by  the  mass  moment 
of  inertia.    This  yields  the  acceleration  in  the  horizontal  and  vertical 
directions  ,  and  the  angular  acceleration  about  the  Y  axis  .    The  accelera- 
tions are  each  integrated  to  give  the  velocities,  and  these,  when  inte- 
grated, give  the  horizontal  and  vertical  position  and  the  angular 
orientation. 

In  this  first  model,  the  tension  force  is  generated  as  a  function  of 
horizontal  velocity  and  body  angle  of  attack.    The  tension  function  was 
formed  in  the  same  manner  as  the  coefficient  equations,  and  was  based 
on  predicted  values  of  the  tension  in  the  cable.    This  tension  function 
also  be  found  in  Subroutine  Tink  of  the  computer  program. 
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The  velocities  and  positions  are  used  to  generate  the  forces  and 
moments  that  are  summed  in  the  basic  equations.    The  forces  and  moments 
on  a  body  are  relative  to  the  velocity  vector  of  the  medium,  and  not  to  the 
local  horizontal.    Thus,  the  sum  of  the  squares  of  the  horizontal  and 
vertical  velocities  yields  the  relative  velocity  squared,  V2  .    The  arc 
tangent  of  the  vertical  velocity  over  the  horizontal  velocity  is  taken  to 
give  the  angle  of  the  velocity  vector,  VAN,  with  respect  to  the  true 
horizontal.     The  angular  orientation  of  the  body  (angle  of  attack),  Al ,  is 
summed  with  VAN  to  produce  the  angle  of  attack  of  the  body  with  respect 
to  the  velocity  vector,  AR;  AR  is  in  radians  .    The  fin  angle  of  attack,  FA, 
is  also  summed  with  VAN  to  yield  the  fin  angle  of  attack  relative  to  the 
velocity  vector,  BR.    These  two  angles  of  attack,  AR  and  BR,  are 
multiplied  by  a  conversion  factor  of  57.29578  deg/rad  to  yield  the  angles 
of  attack,  A  and  B,  in  degrees. 

The  two  angles,  A  and  B,  are  the  arguments  which  determine  the 
coefficient  values  of  body  lift,  drag,  and  moment.     The  coefficients  of 
lift  and  drag  are  both  multiplied  times  V2  and  a  constant  to  produce  the 
body  lift,   BL,  and  the  body  drag,  BD.    The  constant  is  -gpS,  and  is 
equal  to  14.17875. 

The  moment  coefficient  is  multiplied  times  V2  and  a  constant,  f  pSC; 
this  constant  is  equal  to  84.0725.    This  body  moment,  YCM,  is  summed 
with  the  moments  produced  by  the  antenna,  to  yield  the  total  moment, 
acting  on  the  body  about  the  towpoint. 
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There  is  no  definite  information  on  the  antenna,  but  it  is  assumed 
that  it  will  only  contribute  a  drag  force;  this  drag  force  will  contribute  to 
the  total  drag  and  moment  on  the  body.    The  antenna  coefficient  of  drag  is 
assumed  to  remain  constant  at  a  value  of  0.50.    Thus,  the  antenna  drag 
is  dependent  on  a  constant  times  the  frontal  surface  area  times  the 
relative  velocity  squared.    Assuming  a  frontal  width  of  2  inches,  the 
drag  force  of  the  antenna  will  then  vary  with  the  length  of  antenna  under 
water,  and  V2  .    The  length  of  antenna  will  be  approximated  as  the  depth 
of  the  body  below  the  surface.     The  constant  in  this  case  is  ^(l/6)pC       , 
which  equals  0.04146.    The  antenna  drag  is  multiplied  by  the  sine  and 
cosine  of  VAN,  SV,  and  CV,  to  give  two  moment    components.     These  times 
their  respective  moment  arms  yield  YADS  and  YADC . 

The  body  lift  and  drag,  and  antenna  drag  are  multiplied  by  SV  and  CV, 
and  the  resulting  six  values  are  summed  into  the  horizontal  and  vertical 
forces.    In  this  first  realization  of  the  model,  these  horizontal  and 
vertical  forces  are  used  before  being  summed  with  the  horizontal  and 
vertical  components  of  the  tension  force,  to  determine  the  angle  FE  with 
which  the  tension  force  acts  on  the  body.    This  means    that  the  tension 
force  acts  exactly  opposite  to  the  vector  of  body  forces;  this  simplifying 
assumption  is  used  in  the  initial  stages  of  the  model  development. 

The  depth  error  and  depth  rate  inputs  into  the  fin  servo  loop  are  taken 
from  the  vertical  position  and  velocity,   Zand  Zl ,  respectively.     These 
are  multiplied  by  a  constant,   0.443,  which  converts  ft  and  ft/sec  into 
psi/ft  and  psi  sec/ft.    These  are  multiplied  by  gain  constants  ADE  and 
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ADR  before  being  summed  with  the  body  pitch  rate,  Yl ,  times  its  gain 
constant,  APR.     These  three  inputs  are  fed  into  the  fin  servo  loop.    The 
orientation  of  the  fin  causes  a  torque  on  the  body  which  either  increases 
or  decreases  the  body  angle  of  attack;  this  enables  the  body  to  rise  or 
descend. 

This  is  the  initial  dynamic  model;  it  changes  as  more  insight  and 
knowledge  are  gained  on  the  model.     The  first  study  with  the  model 
deals  with  the  three  coefficients  ADE,  ADR,  and  APR. 
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IV.    STUDY  OF  THE  BUOY  MODEL 
Once  the  model  was  digitally  programmed  and  the  initial  programming 
mistakes  ironed  out,  the  first  step  was  to  study  possible  values  of  the  fin 
servo  input  coefficients.    The  first  computer  program  realization  is  shown 
in  the  program  titled  "INITIAL  BUOY  MODEL".    The  buoy  is  initiated  at  a 
depth,   ZI,  of  4.0  ft,  a  speed,  XII  of  20.2536  ft/sec  (12  knots),  body 
angle  of  attack,  All,  of  0.10472  rad  (6  degrees),  and  a  fin  angle  of  attack, 
FAI,  of  0.02  rad  (1.15  degrees).    The  tension  force  and  the  coefficients 
of  lift,  drag,  and  moment  are  called  out  of  Subroutine  Tink,  which  is  in 
double-precision  format. 

A.      SERVO  INPUT  COEFFICIENTS  STUDY 

The  first  attempt  in  this  study  was  to  see  if  the  buoy  would  stabilize 
with  only  one  input  into  the  fin  servo  loop;  either  the  depth  error,   the 
depth  rate,  or  the  body  pitch  rate.     This  proved  fruitless,  so  different 
combinations  of  ADE  and  ADR  were  tried,  and  it  was  found  that  these  two, 
with  APR  set  equal  to  zero,  would  stabilize  the  system.  A  full  range  of 
values  of  ADE  and  ADR  were  tried  to  see  which  combinations  would  bring 
the  buoy  to  a  steady-state  position.     A  mapping  of  these  values  is  shown 
in  Table  IV;  the  value  given  for  each  combination  is  the  depth  at  which 
the  buoy  steadied  out.     The  responses  for  several  of  these  combinations 
are  shown  in  Figs.    12-16.    It  appears  from  the  depth  values  that  the 
tension  force  is  too  strong,  since  in  every  case  the  buoy  steadies  out 
deeper  than  the  desired  depth  of  4.0  ft.    This  became  obvious  when  the 
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steady-state  body  angle  of  attack  and  fin  angle  of  attack  were  checked; 
both  angles  were  large  and  showed  that  the  model  was  working  properly, 
but  that  the  buoy  was  fighting  against  an  unrealistic  tension  force 
formulation. 

It  was  decided  from  the  mapping  of  the  two  coefficients,  ADE  and 
ADR,  that  three  sets  would  be  used  in  further  studies  of  the  buoy  model. 
The  first  pair  was  chosen  on  the  basis  that  it  was  on  the  intersection  of 
the  values  of  ADE  and  ADR  which  had  the  widest  range  of  stable  values, 
and  that  this  set  was  towards  the  center  of  the  mapping.    This  was 
ADE- 1.2  and  ADR=0.30.     It  was  then  decided  that  the  other  two  pairs 
would  vary  in  one  coefficient  away  from  this  point.     These  other  two 
pairs  are  ADE=1  .2  and  ADR=0.07,  and  ADE=2  .5  and  ADR=0  .  30  .    The 
responses  with  these  three  pairs  of  coefficients  are  shown  in  Figs. 
14-16. 

The  coefficient  APR  remained  to  be  studied.    When  this  coefficient 
was  given  a  value  larger  than  zero  there  was  no  appreciable  change  in 
the  response  of  the  buoy.    To  check  the  range  of  possible  values  of  APR, 
the  tension  force,  T,    was  changed  so  that  it  was  equal  to  the  square 
root  of  the  sum  of  the  squares  of  Fl  and  F2  .    Thus,  the  buoy  would  stay 
at  four  feet  below  the  surface     and  the  effect  of  APR  could  be  seen  in  the 
response  of  the  buoy's  angle  of  attack.     Both  negative  and  positive  values 
of  APR  were  tried,  and  negative  values  immediately  drove  the  buoy  moi 
unstable.     Positive  values  of  APR,   up  to  7.0,  were  tried  with  the  three 
pairs  of  depth  error  and  depth  rate  coefficients  .    At  this  upper  value  of 
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7.0,  the  body  started  to  oscillate  away  from  stability  about  the  Y  axis; 
it  was  apparent  that  this  high  value  of  gain  was  approaching  the  stability 
limit  of  the  buoy  model.     Based  on  this,  a  median  value  of  3 . 0  was 
chosen    for  APR.    The  three  final  sets  of  coefficients  are  as  follows: 

ADE  ADR  APR 

1.2  0.07  3.0 

1.2  0.30  3.0 

2.5  0.30  3.0 

B.      TENSION  FORCE  FORMULATION 

The  next  step  in  this  model  study  was  to  formulate  a  new,  and  hope- 
fully more  realistic,  tension  force,  rather  than  using  the  predicted 
tension  data  supplied  with  the  buoy.    It  was  decided  that  this  tension 
function  would  be  a  constant  plus  or  minus  the  distance  away  from  the 
equilibrium  point  times  what  will  be  called  a  spring  constant,  SC  .    The 
basis  for  this  formulation  is  shown  in  Fig.   17.     First  it  is  assumed  that 
the  last  portion  of  the  cable  connected  to  the  buoy,  LO,  is  straight,  and 
that  of  the  total  length  of  the  cable  this  section  is  the  most  dynamic.     By 
computing  the  increase  or  decrease  in   ZO  and  XO,  a  new  cable  length, 
LN,  may  be  calculated.     If  LN  is  larger  than  LO  then  the  tension  increases; 
if  it  is  smaller  then  the  tension  decreases.    The  difference  in  length 
(LN-LO)  is  L  and  this  is  multiplied  times  the  spring  constant,  SC  ,  which 
is  studied  for  several  different  values  later  on.    The  new  tension  force 
now  takes  the  form  : 

T  =  1940  +  L-SC       . 
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THE  NEW  TENSION  FORCE  FORMULATION 
Figure  17 
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The  value  1940  is  the  steady-state  tension  value  that  the  model 
settled  to  in  the  earlier  parameter  studies;  the  angle  FSS  is  the  steady- 
state  value  of  FE.    This  new  tension  force  was  implemented  into  the  buoy 
model.     The  buoy  will  now  be  initiated  at  a  depth  of  zero  feet;  this  is 
possible  because  no  heaving  force  or  surface  effect  is  taken  into  account 
in  the  model;  other  minor  changes  are  made  to  take  this  new  initial 
condition  into  account.     The  new  computerized  model  program  is  titled 
"STUDY  OF  T  =  1940  +  L*SC  " .     The  body  angle  is  now  being  initiated  at 
4.8  degrees  and  the  initial  fin  angle  is  0.0  degrees;  these  two  values 
are  also  the  steady-state  values  for  these  angles  found  earlier.     The  next 
step  in  the  formulation  of  this  tension  force  was  to  find  an  appropriate 
value  for  the  spring  constant.     The  program  was  run  with  the  three  sets 
of  coefficients  and  the  new  initial  conditions  ,  while  SC  was  varied  from 
2000.  to  0.0  in  steps  of  250.     The  two  sets  of  coefficients  with  ADE=1.2 
were  stable  for  all  values  down  to  SC=250,  while  the  set  with  ADE=2.5 
was  stable  down  to  SOI  750.     This  third  set  (ADE=2  . 5 ,  ADR=0.30/  and 
APR=3 . 0)  was  dropped  from  further  investigation  since  this  was  believed 
to  be  an  unrealistic  value  for  the  spring  constant. 

C.     STUDY  OF  INITAL  CONDITION  PERTURBATIONS 

This  final  form  of  the  model  with  the  new  tension  force  formulation 
was  next  examined  to  see  how  stable  the  model  was  under  various  initial 
condition  perturbations  .     It  was  decided  to  perturb  the  initial  vertical 
velocity,  the  body  angle  of  attack,  and  the  fin  angle  of  attack.     Three 
different  spring  constants  were  used  with  the  two  lining  sets  of 
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coefficients;  these  were  SC=500,  SC=1000,  and  SC  =  1500.    It  was  also 
decided  that  two  different  values  for  the  approximation  of  the  dynamic 
length  of  the  cable  would  be  used;  these  were  ZO=100  and  ZO=50.     (ZO 
simply  because  it  was  used  to  initiate  the  model  rather  than  LO) .    The 
final  results  were  tabulated,  and  are  shown  in  Table  V;  the  initial  angles 
of  attack  (All  and  FAI)  are  in  degrees,  and  the  initial  vertical  velocity 
(Z1I)  is  in  ft/sec. 

The  responses  of  the  buoy  model  to  some  of  these  initial  condition 
perturbations  are  shown  in  Appendix  A.     The  twenty-four  plots  include 
all  the  maximum  perturbations  listed  for  SC=500  and  SC  =  1500.    These 
plots  do  not  really  tell  much  about  the  buoy  model.     They  do  show  that 
the  buoy  will  settle  to  steady  state  fairly  rapidly,  and  that  the  buoy  will 
not  deviate  very  far  from  its  desired  depth.     Plots  with  the  same  pair  of 
coefficients  (ADE  and  ADR),  with  the  same  dynamic  length  of  cable,  and 
with  the  same  type  of  inital  condition  perturbation,  but  with  different 
values  of  spring  constant,  seem  to  be  very  similar.     The  responses  with 
the  larger  spring  constant  (SC=1500)  seem  to  settle  more  quickly,  even 
though  they  are  oscillating  through  a  wider  range  of  depths  .    All  the 
responses  to  the  fin  angle  perturbations  are  remarkably  similar.     This 
might  be  because  the  fin  is  quickly  driven  back  to  the  desired  stability 
position;  and  since  all  the  fin  angle  perturbations  are  the  same  pair  of 
values,  this  similarity  in  responses  is  not  surprising.    The  perturbing  of 
the  initial  conditions  for  the  new  tension  force  realization  was  the  [ 
at  which  this  study  ended. 
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ZO=50.  ZO=100. 

ADE=1.2  ADE=1.2  ADE=1.2  ADE=1.2 

ADR=0.07  ADR=0.3  ADR=0.07  ADR=3 . 0 

APR=3.  APR=3.  APR=3.  APR=3 . 

SC=500 

All  -2.1-+0.9  -1.8-+0.9  -1.7-+1.  -1.5-+0.9 

All  +2.5-+8.0  +3.0-  +7. 5  +2.0- +7.  +2. 5- +6. 5 

FAI  -17.5-+15.0  -17.5-+15.0  -17.5-+15.0  -17.5-+15.0 


SC=1000 

Z1I 

-4.8-+2.8 

-2.9-+2.6 

-4.3-+2.5 

-2.8-+2.2 

All 

-0.5-+9.5 

0,0-  +9.0 

-0.5-+8.0 

0.0-+8.0 

FAI 

-17.5-15.0 

-17.5-+15.0 

-17.5-+15.0 

-17.5-+15. 

SC  =  1500 

Z1I 

-4. 8- +3.1 

-2.7-+3.3 

-4.8-+3.0 

-2.7-+2.9 

All 

-2.0-+10.0 

-1  .0-+9.5 

-1.5-  +9.0 

-1.0-+8.5 

FAI  -17.5-+15.0  -17.5-+15.0        -17.5-+15.0       -17.5-+15.0 

STABILITY  RANGE  OF  INITIAL  CONDITIONS 
Table  V 
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V.     CONCLUSIONS  AND  RECOMMENDATIONS 
The  study  presented  here  is  only  a  beginning  to  this  problem.     There 
are  many  facets  and  additions  to  this  problem  that  remain  to  be  studied. 
First  of  all,  the  dynamic  characteristics  and  wave-following  capabilities 
of  this  model  were  never  investigated.     Other  additions  could  be  made  to 
the  model  to  take  into  account  the  heaving  forces  generated  by  the  waves 
and  the  free-surface  effect  which  will  definitely  affect  the  motion  and 
capabilities  of  the  buoy.    There  remains  much  to  be  studied  in  the  area 
of  the  dynamic  cable  and  the  coupling  of  the  tension  force  to  the  buoy 
model.     Other  tension  force  realizations  might  be  tried  with  several  spring 
constants  in  a  single  cable  length;  or  the  cable  might  be  approximated  by 
a  series  of  straight  lengths,  each  with  its  own  spring  constant.    The 
possibilities  of  adding  a  damping  coefficient  into  the  tension  force  real- 
ization need  to  be  investigated;  the  possibilities  are  almost  endles  . 
Other  methods  of  simulating  the  buoy  need  to  be  tried,  especially  that 
of  using  a  hybrid  computer  (an  analog  computer  integrated  with  a  digital 
computer);  hybrid  computers  are  more  suited  to  this  type  of  real-time 
simulation  work.    A  dynamic  model  that  allowed  for  all  the  possible 
dynamic  forces  and  all  the  possible  variations  in  buoys  and  cables, 
would  be  invaluable  in  the  design  and  engineering  of  towed  buoys  . 


APPENDIX  A 
RESPONSES  OF  THE  BUOY  TO  INITIAL  CONDITION  PERTURBATIONS 
The  following  24  plots  are  the  responses  of  the  buoy  to  the  maximum 
initial  condition  perturbations  listed  in  Table  V  for  SC=500  and  SC=1500. 
The  plots  are  of  depth  versus  time.    The  x  scale  for  all  the  plots  is 
5  sec/inch;  and  the  y  scales  are  in  ft/inch. 
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APPENDIX  B 

SYMBOLS  AND  ABBREVIATIONS  USED  IN  THE  COMPUTER  PROGRAMS 

Z  vertical  position  in  ft.   (initial  condition  ZI) 

Zl  vertical  velocity  in  ft/sec.    (initial  condition  Z1I) 

Z2  vertical  acceleration  in  ft/sec2 

FV  summation  of  vertical  forces 

Fl  summation  of  the  vertical  forces  less  the  vertical  tension 

component 

XI  horizontal  velocity  in  ft/sec.   (initial  condition  XII) 

X2  horizontal  acceleration  in  ft/sec2 

FH  summation  of  the  horizontal  forces 

F2  summation  of  the  horizontal  forces  less  the  horizontal 

tension  force 

FF  the  absolute  value  of  Fl 

FG  the  absolute  value  of  F2 

FE  the  angle  with  which  the  tension  force  acts  on  the  body 

in  radians 

SFE  the  sine  of  FE 

CFE  the  cosine  of  FE 

TSF  the  vertical  component  of  the  tension  force 

TCF  the  horizontal  component  of  the  tension  force 

C  horizontal  velocity  in  knots 

VAN  the  angle  of  the  velocity  vector  relative  to  the  horizontal 

in  radians 

V2  the  magnitude  of  the  velocity  vector  squared 
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Al  the  body  angle  of  attack  relative  to  the  horizontal  in  radians 

(initial  condition  All) 

Yl  the  body  pitch  rate  about  the  Y  axis  in  radians/sec 

(initial  condition  Yll) 

Y2  angular  acceleration  of  the  body  about  the  Y  axis  in 

radians/second2 

YM  summation  of  the  moments  on  the  body  about  the  Y  axis 

. 

AR  the  body  angle  of  attack  relative  to  the  velocity  vector 

in  radians 

A  AR  in  degrees 

FA  the  fin  angle  of  attack  relative  to  the  horizontal  in 

radians  (initial  condition  FAI) 

FB  used  in  the  fin  servo  loop 

FC  used  in  the  fin  servo  loop 

FD  summation  of  the  inputs  and  the  feedback  in  the  fin 

servo  loop 

BR  the  fin  angle  of  attack  relative  to  the  velocity  vector  in 

radians 

B  BR  in  degrees 

CL  coefficient  of  lift  of  the  body 

CD  coefficient  of  drag  of  the  body 

CM  coefficient  of  moment  of  the  body 

T  the  tension  force 

CV  the  cosine  of  VAN 

SV  the  sine  of  VAN 

ZAB  the  absolute  value  of  z 

AD  the  antenna  drag 
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BL  the  body  lift 

BD  the  body  drag 

YCM  the  moment  on  the  body  about  the  Y  axis 

YADS         the  moment  on  the  body  due  to  the  vertical  component  of 
the  antenna  drag 

YADC         the  moment  on  the  body  due  to  the  horizontal  component 
of  the  antenna  drag 

ZD  in  the  first  program  this  is  used  to  measure  the  relative 

depth  error  from  4.0  ft;  in  the  second  program  it  is  used  to 
compute  the  antenna  drag 

DEE  the  depth  error  in  psi/ft 

DE  DEE  times  the  coefficient  ADE 

DRR  the  depth  rate  in  psi- sec/ft 

DR  DRR  times  the  coefficient  ADR 

YPR  the  body  pitch  rate  (Yl)  times  the  coefficient  APR 

ADE  the  depth  error  coefficient 

ADR  the  depth  rate  coefficient 

APR  the  pitch  rate  coefficient 

THE  FOLLOWING  ARE  USED  SOLELY  IN  THE  SECOND  COMPUTER  PROGRAM 
IN  THE  REALIZATION  OF  THE  NEW  TENSION  FORCE 

FSS  the  steady-state  value  of  FE 

ZO  the  vertical  approximation  to  the  dynamic  length  of  the 

cable;  see  Fig.   17 

SC  the  spring-constant  value  of  the  cable  in  pounds/ft 

SFS  the  sine  of  FSS 

LO  the  approximation  to  the  dynamic  length  of  t)  Le; 

see  Fig .   1 7 


TFS  the  tangent  of  FSS 

XO  the  horizontal  approximation  to  the  dynamic  length  of  the 

cable;  see  Fig .   17 

ZN  the  new  computed  value  of  ZO 

XR  the  relative  movement  out  of  steady  state  in  the  X 

direction 

X  the  change  in  position  out  of  the  steady  state  in  the 

X  direction 

XN  the  new  computed  value  of  XO 

LN  the  new  computed  length  of  the  dynamic  length  of  the 

cable 

L  the  change  in  length  of  the  cable 

FE  as  stated  previously,  but  computed  differently 
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